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The effect of the 19 F(a, p) 22 Ne reaction rate uncertainty on 
the yield of fluorine from Wolf-Rayet Stars 
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ABSTRACT 

In the light of recent recalculations of the 19 F(a,p) 22 Ne reaction rate we present 
results of the expected yield of 19 F from Wolf-Rayet (WR) stars. We have computed 
models using the upper and lower limits for the rate in addition to the recommended 
rate and hence we constrain the uncertainty in the yield with respect to this reaction. 
We find a yield of 3.1 x 1CP 4 M0 of 19 F with our recommended rate and a difference 
of a factor of two between the yields computed with the upper and lower limits. In 
comparison with previous work we find a difference in the yield of approximately a 
factor of 4, connected to a different choice of mass loss. Model uncertainties must be 
carefully evaluated in order to obtain a reliable estimate of the yield of fluorine from 
WR stars together with its uncertainties. 

Key words: stars: evolution, stars: Wolf-Rayet, stars: interiors, nucleosynthesis, 
fluorine 



1 INTRODUCTION 

The production site of 19 F has been a major p uzzle for nucle- 
osynth esis for a long time. It was predicted bv lGorielv et alJ 
( 1989) that 19 F should be manufactured in asymptotic giant 
branch stars, and these are currently the onl y observation- 
ally confirmed site for fluorine production jjorissen et alJ 
1992). Other sites and mechanisms for the galactic produc- 
tion of fluorine have been proposed. The neutrino process 
operatin g during type-II supernov a explosions can produce 
fluorine iWooslev fc Weaver! Il995h . Moreover, fluorine can 
be synthesised during core He-burning and ejected via the 
strong wind s of Wolf-Rayet (WR) stars dMevnet fc Arnouldl 
Il99l hood) . It appears that the contributions of asymp- 
totic giant branch and WR stars have to be included in 
the computation of the chemical evolution of the Galaxy to 
account for the ob servations of fluorine in the Milky Way 
l|Renda et alJl2004) . 

The uncertainties of the nuclear reaction rates involved 
in the production and destruction of 19 F in the case of 
as ymptotic giant branc h stars have been studied in detail 
bv lLuearo etal] (I2004D . Here we examine fluorine produc- 



E-mail: rs@ast.cam.ac.uk 



tion in WR stars with regard to the uncertainties of the 
19 F(a,p) 22 Ne rate. 

Fluorine production in WR stars is of secondary nature 
because it relies on the presence of 13 C and 14 N produced 
during H burning. During core helium burning the produc- 
tion mechanism for 19 F is as follows. The seed 14 N under- 
goes a capture to produce the unstable nucleus 18 F which 
/3-decays to 18 O. In the presence of protons the reaction 
ls O(p, a) 15 N can occur which then leads to the formation of 
19 F via the reaction 15 N(q, 7) 19 F. Protons may be present in 
a helium burning core from the reaction 14 N(n, p) 14 C, with 
the neutrons coming from the reaction 13 C(a, n) ie O. The 
important destruction mechanism associated with the core 
helium burning phase is the reaction 19 F(cv, p) 22 Ne. New es- 
timates for the rate of this reaction have been presented by 
iLugaro et alJ (2004) . It is very difficult to gather experimen- 
tal data at low energies so the estimated uncertainties are 
of fourteen orders of magnitudes at a temperature of about 
2 x 10 s K, typical of core helium burning. Thus it is neces- 
sary to investigate the impact of such a large uncertainty on 
the yield of fluorine from WR stars. 

In section|2]we describe the details of the evolution code 
used and its nucleosynthesis routine. We briefly review the 
details of the 19 F(tv, p) 22 Ne reaction in sectional In section^] 
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Light isotopes 



Reaction 



Source 



1 H,n 1 2 H, 3 He, 4 He 7 Li, 7 Be, 11 B, 12 C, 13 C, 
^C, 14 N, 15 N, 16 O, 17 O, 18 O, 19 F, 20 Ne, 21 Ne, 22 Ne, 
8«Aro, 23 Na, 24 Mg, 25 Mg, 26 Mg, 2fl Ai m , S6 ^9, 27 Al, 28 Si, 29 Si, 

30gj 31p 32g 33g 34g 



Iron group isotopes 



56 Fe, 57 Fe, 58 Fe, 59 Fe , 60 Fe 
59 Co, 58 Ni, 59 M, 60 Ni, 61 Ni 



Table 1. Isotopes included in the nucleosynthesis code. Isotopes 
also included in the structural part of the code are highlighted in 
bold. Unstable isotopes are in italics. 



we present the results of the simulations. Conclusions and 
directions for further work are outlined in section |5] 



2 THE STELLAR EVOLUTION CODE 

Calculations were made using the evol ution code STARS 
originally developed by lEggletonl Jl97lf) and updated by 
many contributors (e.g. |Polset^j]|l995l) . The current ver- 
sion of the code employs a fully simultaneous solution of 
the equations of stellar structure, nuclear burning and mix- 
ing ([S tancliffe. Tou t fe Polsll2004fl . Opacities are taken from 
iRogers fe Iglesiasl il992T ) and lAlexander fe Ferguson! lll994) . 
The evolution code uses an approximate reaction network 
comprising only those elements - namely 1 H, 4 He, 12 C, 14 N, 
le O and 20 Ne - whose reactions are energetic enough to af- 
fect the stellar structure. 

In order to fully investigate the nucleosynthesis that 
takes place in stars we require an additional routine that 
comprises a full network for the elements that we are in- 
terested in. Within the framework of the STARS code, this 
takes the form of a subroutine called after the convergence of 
a model with the main code. The nucleosynthesis subroutine 
takes the structural details from this converged model and 
uses them to compute the nucleosynthesis of all the elements 
we are interested in. 



2.1 The nucleosynthesis routine 

The nucleosynthesis routine deals with 40 isotopes. These in- 
clude all stable particles, and a few unstable ones, from neu- 
trons and deuterium up to 34 S and some isotopes from the 
iron group. It also computes its own values for the isotopes 
included in the evolutionary code as the compositions from 
the evolutionary part of the code will deviate slightly from 
those of the nucleosynthesis part. A list of all the isotopes 
included in the nucleosynthesis code is given in Tabled 

Unstable nuclei that are not included in the network 
are treated as if their decay were instantaneous. This ap- 
proximation is fair for all light isotopes with half-lives from 
seconds to hours in conditions of core helium burning. For 
the unstable isotopes considered in the ne twork, t he de cay 
lifetimes are the terrestrial values given bv iKrana il988f) . 
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Table 2. Proton capture reactions a nd the sources from whic h 
their rates were taken. Key: C88 JChamgagne^^alJ ^983) , 
CF88 iCaughlan fe Fowled Il988t) . B95 jBlackmon et al.l Il995l) . 
EL95 <EI Eid fc Champagnelll995l). G90 Icprres et al.ll990l). 190 
illiadis et alJll990j). 196 llliadis et alj|l99rt. L90 jLandre et alj 
I1990J) , NACRE lAngulo et alll999l) . SC95 jschmidt et al.ll995l) . 
S T96 iStegmuller et alJll996t) , REACLIB (1991 updated v ersion 
oflThielemann et al.ll98rj) . T88 iTimmermann et al.ll98sT) . W90 
iWiescher et alJll99oE 



2.1.1 Charged particle reaction rates 

In order to couple the nucleosynthesis network 63 charged 
particle reactions are required. The rates are taken from a 
variety of sources and are listed in Table (proton cap- 
tures) and Table [3] (a captures ). The rate of the reactio n 
3 He( 3 He,2p) 4 He is that given bv lCaughlan fe Fowled Jl98d) . 
as are the rates for carbon and oxygen burning reactions. 

The nucleosynthesis routines were designed to employ 
the ready-to-use fits to the reaction rat es from the REA- 
CLIB library (1991 updated version of iThielemann et alJ 
1986), updated wher e possible to includ e the latest exper- 
imental results (see iLugaro et alJ 120041 for full details) . 
For some of the rates involved in the production of 19 F, 
such as 15 N(a,7) 19 F, the rates are virtually the same as 
those presented in the NACRE compilation jAngulo et alJ 



1999). In other cases, such as the rates N(a,7) F and 
18 0(a, 7) 22 Ne, the rates used are updates with respect to 
the NACRE rates. 
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Reaction 


Source 


4 He(aa, 7 ) 12 C 


CF88 


7 Li(o,7)llB 
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12 C(a,7) 16 
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14 C(a, 7 ) 18 
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14 N(a,7) 18 F 


GOO 


15 N(q, 7 ) 19 F 


de096 


16 O(a,7) 20 Ne 


CF88 


17 O(a,n) 20 Ne 


D95 


18 0(a,7) 22 Ne 


D03 


18 0(a,n) 21 Ne 


D95 


19 F(a,p) 22 Ne 


U04 


20 Ne(a,7) 24 Mg 


CF88 


21 Ne(a,7) 25 Mg 


CF88 


21 Ne(a,n) 24 Mg 


D95 


22 Ne(a, 7 ) 26 Mg 


K94 


22 Ne(a,n) 25 Mg 


K94 


^Nafa, n) 2b A19 


CF88 


23 Na(a,n) 26 Al m 


CF88 


24 Mg(o, 7 ) 28 Si 


CF88 


25 Mg(a, 7 ) 29 Si 


CF88 


25 Mg(a,n) 28 Si 


CF88 


25 Mg(«,p) 28 Al 


CF88 


26 Mg(a,7) 30 Si 


CF88 


2e Mg(a,n) 29 Si 


CF88 


26 Mg(a,p) 29 Al 


CF88 


27 Al(a, 7 ) 31 Si 


REACLIB 



Table 3. Reactions involving a capture a nd the sources from 
which their rates were taken. Key: CF88_j2S!JSiii a J^^iffiiSl 
Il988l). D95 toenker et alj|l995l) . deQ96 <de Oliveira et al.ll996l) . 
GOO iGorres et al\\2()Kt . ,TG 01 i.Torissen Gorielvll200ll) . D03 
iDababneh et all2003ft . K94 iKaeppeler et al.ll994l) . RE ACLIB 
1991 updated version of iThielemann et al J <198fil) . U04 iUgaldel 
l2004h . 



2.1.2 Neutron capture rates 

A total of 45 neutron capture reactions a re req uired 
for the network. The work of iBao et all i2000t) was 
used as th e main source. Supplem e ntary (n, 7) data are 
taken from iRauscher fc Thielemann! ()200(tl for captures by 
59 ' 60 F e . Rates for the reaction 33 S(n, a) 30 Si were taken 
from lSchatz et alJ lll995l) . For the rea ctions 26 Al 9 (n p) 26 M g 
and 26 Al 9 (n,a) 23 Na rates are from iKoehler et alJ il997f) . 
The important react ion rate for 14 N(n, p) 14 C is from 
Glcd enov et al ] lll995l). whic h is in agreeme nt with previ- 
ous experimental l|K oehler fc 0'brienlll989T) and theoreti- 
cal (iBahcall fc Fowlerlll969f estimates. This rate is approx- 
imately a factor of two higher t han the rate propo s ed by 
[ Brehm et all Jl988l) and used bv lMevnet fc Arnouldl J19951 

For neutron capture s by 59 Ni we take reaction rates 
from lHolmes et ail lll976l) and this is also the source of the 
rate of the reaction lr O(n, a) 14 N. 

In addition to this, two neutron sinks are included to 
account for neutron captures by those elements not included 
in the network. The first sink is emulated by the reac- 
tion 34 S(n, 7) 35 S and represents nuclei between 34 S and the 
iron group. The second sink is emulated by the reaction 
61 Ni(n, 7) 62 Ni and represents captures by all the heavy ele- 
ments above 61 Ni. In any case these reactions are not very 
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Figure 1. The recommended rate of the 19 F(o, p) 22 Ne reaction 
is plotted as a solid line. Upper and lower limits of the rates 
are plotted as dotted lines, while the dashed line represents the 
22 Ne(a, n) 25 Mg reaction rate. 

important in core He-burning conditions where the main 
neutron absorber is 14 N. 



3 THE 19 F(a,p) 22 Ne REACTION RATE 

The 19 F(a,p) 22 Ne reaction rate is the main destruction 
channel of 19 F in WR stars. However there is very little 
experime ntal data avai l able a t low energies. As described in 
detail bv lLugaro et alJ (|2004lh . the level density in the com- 
pound nucleus 23 Na has been analysed on the basis of, as yet 
un published, low- energy 19 F(a,p) resonance measurements 
bv llJgaldel l)2004h . We found that the leve l density is too 
low t o apply the Hauser-Feshbach approach iRauscher et alJ 
119971) . whic h yields a rate in reasonabl e agreement with the 
estimate bv lCauehlan fc Fowled (I1988T) . Thus the rate needs 
to be calculated from determination of the strengths, u>"/, 
for the single resonances. The resulting recommended rate 
is shown in Figure together with the upper and lower lim- 
its. In the temperature range of core He burning the recom- 
mended rate is more than one order of magnit u de sm aller 
than the rate estimated bv lCauehlan fc Fowled lll988t) and 
used in the previous calculations. The lower limit for the 
rate is several orders of magnitude smaller than the recom- 
mended rate. 



4 RESULTS 

To test the effects of varying the 19 F(a, p) 22 Ne reac- 
tion rate a 60 Mq model of Z = 0.02 was evolved 
from the pre-main sequence through to the WR phase 
with the recommended reaction rate and the upper and 
lower limits. Initial abu ndances for all the e lemen ts con- 
sidered are taken from lAnders fc Grevessd (I1989T) . This 
means we take the initial 19 F mass fraction as 4.051 x 



4 



R.J. Stancliffe, M. Lugaro, C. Ugalde et al. 



Reaction Rate Yield (10 -4 M©) 



5.5e+08 
5e+08 
4.5e+08 
£ 4e+08 
| 3.5e+08 
2 3e+08 
| 2.5e+08 
u 2e+08 
° 1.5e+08 
le+08 
5e+07 


4e+06 4.1e+06 4.2e+06 4.3e+06 4.4e+06 4.5e+06 4.6e+06 
Time (years) 

Figure 2. Core temperature as a fuction of time for the He- 
burning phase. 



10~ 7 . Mass loss was applied by the ma s s-loss rates of 
Ide Jager. Nieuwenhuiizen fc van der Hu chd (Il98fif) for pre- 
WR evolution and the rates of lLangerl IU989T) for the WR 
phase. The switch to WR mass loss is made when the surface 
hydrogen abundance by mass fraction reaches 0.4 and the 
surface temperature exceeds 10 4 K. These may no t be the 
ideal prescriptions to use fsee lEldridee fc Toutfeool but our 
aim is to investigate the effect of varying the 19 F(a, p) 22 Ne 
reaction rate not the influence of mass-loss rate on WR evo- 
lution. 

Figure |H] shows the mass fraction of 19 F in the core 
and at the surface as a function of time. Initially a rapid 
increase in the core fluorine abundance occurs because all 
the available 13 C and 14 N undergo a captures, opening the 
path to 19 F production. It takes about 8 x 10 4 yr from this 
time to have a notable increase in the abundance of 19 F at 
the surface when stellar winds have stripped away the stellar 
surface exposing the formerly convective core. Between this 
point and the end of the life of the star (the WC phase) 
approximately 9M© of material is lost from the surface. 

A plot of the core temperature as a function of time dur- 
ing the He-burning is shown in Figure [5] In the case of the 
upper limit of the 19 F(a, p) 22 Ne reaction, we see that tem- 
peratures are sufficiently high for the reaction to be active in 
the core soon after the maximum abundance of 19 F has been 
reached and the fluorine abundance in the core rapidly drops 
from its peak value. For the recommended rate the initial 
decline is much slower and becomes steeper after 4.4 x 10 6 
yr. The lower limit to the rate leads to an almost constant 
core abundance of 19 F until about 4.4 x 10 6 yr later at which 
point fluorine is appreciably destroyed. The destruction of 
fluorine in the lower-limit case, and the steeper decline in the 
recommended case, that occur at a time of about 4.4 x 10 6 yr 
are to be attributed to another destruction channel that is 
opening at around such a time in the centre of the star, the 
19 F(n, 7) reaction where the neutrons are provided by the 
activation of the 22 Ne(a, n) 25 Mg reaction rate. The rate we 
use for this neutron source reaction is practically the same 
as that recommended by NACRE. Such a rate is comparable 
to the recommended limit of the 19 F(a, p) 22 Ne reaction and 
it is much higher than the lower limit for this reaction, as 
illustrated in Figure The estimate of the neutron capture 
cross section of 19 F is relatively old and needs to be revalu- 



Upper limit 
Recommended 
Lower limit 



1.7 

3.1 
3.4 



Table 4. The yield of 19 F obtained with the indicated 
19 F(a, p) 22 Ne reaction rates. 



ated feao et al . 2000) but a different rate is unlikely to have 
a large impact on the 19 F yield from WR stars as the sur- 
face abundance is not affected by this destruction channel, 
as shown by the left panel of Figure [3] 

The yield of 19 F from each of the evolution runs is pre- 
sented in Table 2] It is denned as 



p= I M(t)[X s (t)-X ]dt 



(1) 



where r is the lifetime of the star, M(t) the mass-loss rate 
at age t, X s (t) the surface fluorine abundance at time t 
and Xo the initial abundance of fluorine. The rapid de- 
struction of fluorine that occurs with the upper limit of the 
19 F(«,p) 22 Ne rate leads to a yield approximately a factor 
of two lower than the cases computed using the recommend 
and the lower limit for the rate. 

In order to compare our resu lts with those presented by 
iMevnet fc Arnouldl dl99l l2000l) for the equivalent model 
we must consider the results of our calculation computed 
with the upper limi t for the rate. This is eq uivalent to 
the rate proposed bv lCauehlan fc Fowled 1^988). Our yield 
is very close to the value of 1.2 x 10 -4 M Q calculated by 
IMevnet fc Arnouldl lll99&tl , who used the standard mass-loss 
rate that we adopted. This yield is ap proximately a factor 
of fou r smaller than that calculated by IMevnet fc Arnouldl 
(2000), who adopted an enhanc ed mass-loss rate. The en - 
hanced mass loss rates used by IMevnet fc Arnouldl ll200Cf) 
will expose the He core sooner and so provide a higher 19 F 
yield. The comparison underlines the importance of mass- 
loss rates on the calculations of yields from these stars. 

We also note that many reaction rates that influence 
the production of fluorine have been updated in our calcu- 
lations with respect to those of ICauehlan fc Fowlerl (ll98ST) 
that were used by Meynet fc Arnould. Moreover, we have 
not included any convective overshooting. This means that 
we have a smaller convective core than Meynet fc Arnould, 
which affects the luminosity of the star and thus the mass- 
loss rate. 



5 CONCLUSIONS 

We have calculated the yield of 19 F from a Wolf-Rayet star 
of 60 Mq at Z = 0.02 with the upper, recommended and 
lower limits for the rate of the 19 F(cv, p) 22 Ne reaction. We 
find that there is a difference of a factor of two between 
the yields computed with the upper and lower limits for the 
rate. As future work the effects of the new rates for the 
19 F(«,p) 22 Ne reaction need to be determined across a wide 
range of stellar masses and metallicities and the effect of 
the uncertainties associated with the other rates affecting 
the production of fluorine should also be analysed. 

We also find that, for a given rate of the 19 F(a,p) 22 Ne 
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Figure 3. Left panel: mass fraction of 19 F present in the core as a function of time. The sharp increase in the mass fraction at 4.1 X 10 6 
yr is the onset of core He-burning. Right panel: mass fraction of 19 F at the stellar surface as a function of time. 



reaction, the F from Wolf-Rayet stars is reduced by a fac- 
tor of about 4 with respect to p revious calculations, mad e 
with an enhanced mass-loss rate (iMevnet fc Arnoulcl20o"o|) . 
These model uncertainties should be carefully evaluated. 
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